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HARD X-RAY EMISSION FROM ACCRETION SHOCKS AROUND GALAXY CLUSTERS 
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ABSTRACT 

We show that the hard X-ray (HXR) emission observed from several galaxy clusters is consistent with 
a simple model, in which the nonthermal emission is produced by inverse Compton scattering of cosmic 
microwave background photons by electrons accelerated in cluster accretion shocks: The dependence 
of HXR surface brightness on cluster temperature is consistent with that predicted by the model, 
and the observed HXR luminosity is consistent with the fraction of shock thermal energy deposited 
in relativistic electrons being < 0.1. Alternative models, where the HXR emission is predicted to 
be correlated with the cluster thermal emission, are disfavored by the data. The implications of our 
predictions to future HXR observations (e.g. by NuStar, Simbol-X) and to (space/ground based) 
7-ray observations (e.g. by Fermi, HESS, MAGIC, VERITAS) are discussed. 

Subject headings: acceleration of particles - galaxies: clusters: general - radiation mechanisms: non- 
thermal - X-rays: general 



1. INTRODUCTION 

Nonthermal emission is observed from several clus- 
ters of galaxies, mainly in the radio band (e.g., Feretti 
& Giovannini 2008). In some cases, nonthermal hard 
(> 20keV) X-ray (HXR) emission is also observed (for 
review, see Rephaeli et al. 2008). The radio emission is 
interpreted as synchrotron radiation, thereby suggesting 
that relativistic electrons and magnetic fields are present 
in the intracluster medium (ICM). The HXR emission 
is usually interpreted as due to inverse Compton (IC) 
scattering of cosmic microwave background (CMB) pho- 
tons by nonthermal relativistic electrons (e.g., Rephaeli 
1979; Sarazin 1999). In some cases, however, the HXR 
emission is also consistent with a two-temperature ICM 
plasma (see Rephaeli et al. 2008, and refrences therein). 

An accurate determination of the nonthermal HXR 
flux requires a precise measurement of the thermal emis- 
sion, which dominates up to ~ 30 — 40 keV. This, as well 
as a clear answer to the question of whether the HXR 
emission is due to a two-temperature plasma or to the 
presence of a non-thermal electron population, will only 
be provided by future HXR and 7-ray missions. In this 
paper, we examine the consequences of the nonthermal 
HXR detections reported in the literature, assuming that 
the reported fluxes are correct (i.e. that the substraction 
of the thermal component has been preformed correctly) 
and that the energy distribution of the nonthermal pho- 
tons is well described by a power-law (and is not due, 
therefore, to a two-temperature plasma). We consider 
all reported measurements of HXR emission, excluding 
those where confusion with AGNs or radio galaxies (lo- 
cated within the field of view of the cluster observation) 
is possible (At present, data are available for 13 clusters 
satisfying this criterion, see § 2 for details). 

Several models for the HXR emission from clusters 
have been presented in the literature. These models dif- 
fer in the assumptions regarding the nonthermal emis- 
sion mechanism as well as regarding the origin of the 
emitting electrons. In some models, the nonthermal 
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emission mechanism is IC scattering of CMB photons 
by relativistic electrons (e.g., Rephaeli 1979; Sarazin 
1999; Colafrancesco & Marchegiani 2009), while in oth- 
ers the mechanisms are nonthermal bremsstrahlung (e.g., 
Sarazin 1999; Sarazin & Kempner 2000) or synchrotron 
emission from ultra-relativistic electrons (Timokhin et al. 
2004; Inoue et al. 2005; Eckert et al. 2007b). Various 
sources have been suggested for the emitting electrons: a 
population of point sources (e.g. AGN as in Katz 1976; 
Fabian et al. 1976; Fujita et al. 2007), merger shocks 
(e.g., Fujita et al. 2003; Brunetti et al. 2004), dark mat- 
ter bow shocks (e.g., Bykov et al. 2000), ram-pressure 
stripping of infalling galaxies (e.g., de Plaa et al. 2006) 
and accretion shocks (e.g., Loeb & Waxman 2000; Fujita 
et al. 2003; Berrington & Dermer 2003; Gabici & Blasi 
2003; Brunetti et al. 2004; Inoue et al. 2005; Kushnir & 
Waxman 2009). 

We examine in this paper the consistency of the predic- 
tions of different models for the HXR emission with avail- 
able HXR observations. Two types of models are consid- 
ered. In models of the first type, the HXR emission is 
predicted to be correlated with the cluster thermal emis- 
sion, and therefore to be strongly concentrated towards 
the cluster's center. Models of this type are largely moti- 
vated by the spatial correlation observed in some clusters 
between the nonthermal radio emission and the thermal 
X-ray emission. A widely discussed model of this type is 
a model where the HXR emission is due to secondary 
electrons produced by inelastic p-p collisions between 
cluster cosmic-rays (CRs) and thermal intra-cluster gas 
(e.g., Dennison 1980). In this model, the nonthermal 
radio emission and the HXR emission are produced by 
the same population of electrons. In the second type of 
models, the nonthermal HXR emission originates in the 
cluster accretion shocks, and is therefore extended across 
the cluster. In such models the radio and HXR radiation 
are produced by different electron populations. 

The cluster sample that we use is described in § 2. 
In § 3 we show that all available HXR cluster observa- 
tions are consistent with a simple analytic model (Loeb & 
Waxman 2000; Kushnir & Waxman 2009), in which the 
HXR emission is due to IC scattering of CMB photons 
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by relativistic electrons accelerated in accretion shocks 
surrounding the clusters (§ 3 includes a brief description 
of the model. The reader is referred to Kushnir & Wax- 
man 2009, for a detailed description). In § 4 we show 
that models, in which the HXR emission is predicted to 
be correlated with the cluster thermal emission, are dis- 
favored by the data. This conclusion is based on two 
main findings. First, we show that Swift's upper lim- 
its on HXR emission from several clusters are difficult 
to explain in models where the emission is dominated 
by the cluster core (while naturally explained in the ac- 
cretion model). Second, we show that in models where 
HXR emission is due to secondary electrons produced by 
inelastic p-p collisions, the energy in CR protons is re- 
quired to exceed the thermal energy of the gas in order to 
explain the detected fluxes. This is both unlikely and in- 
consistent with the observed correlation between the ra- 
dio flux and the thermal X-ray flux (Kushnir et al. 2009). 
Our results are summarized and discussed in § 5, with 
emphasis on predictions of the accretion model that dis- 
criminate it from other models and which may be tested 
by future observations. 

The following point should be clarified here. It is diffi- 
cult to directly determine using current data whether the 
observed HXR emission is extended or dominated by the 
cluster core. This is due mainly to the fact that imaging 
information in the HXR band is at best limited 2 . We 
show, however, that this difficulty may be partially over- 
come, and that information regarding the spatial distri- 
bution of the HXR emission may be obtained, by compar- 
ing the fluxes measured by different instruments which 
differ in their fields of view (FOV). 

The FOV radii of RXTE, BeppoSAX and INTEGRAL 
arc in the range of 30' to 60', and those of Swift/BAT 
and Suzaku are smaller, ~ 10'. Since these FOV radii are 
comparable to the characteristic sizes of massive clusters 
(lying at distances of few 100 Mpc), the HXR flux mea- 
surements obtained by instruments with differing FOVs 
provide some information on the intra-cluster spatial dis- 
tribution of the HXR emission. Models in which the 
HXR emission follows the thermal emission predict that 
the HXR flux should not significantly vary as the FOV 
grows to include cluster regions beyond the cluster core, 
since the thermal emission is strongly dominated by the 
core. On the other hand, models in which the HXR sur- 
face brightness is roughly uniform across the cluster, or 
rising away from the cluster center (as predicted by the 
model described in § 3), predict that the HXR flux should 
increase significantly as the FOV grows beyond the angu- 
lar size of the core. As we show in § 4, the data support 
the latter qualitative behavior, therefore suggesting that 
the HXR emission is extended and not dominated by the 
cores of the clusters. It should, however, be kept in mind 
that comparing the fluxes measured by different instru- 
ments is subject to uncertainties, since different instru- 
ments may be subject to different systematic effects (see, 
e.g., Rossetti & Molendi 2004, 2007). Thus, only future 
HXR (and 7-ray) missions, capable of producing high 
resolution maps of clusters, would provide a clear deter- 

2 One exception is the INTEGRAL measurements of HXR emis- 
sion from the Coma cluster, which indicate a source extended well 
beyond the cluster core (Eckert et al. 2007a; Lutovinov et al. 2008). 
A firm conclusion can not, however, be drawn, since the statistical 
significance of the Coma HXR detection by INTEGRAL is low. 



mination of the spatial distribution of the HXR emission 
(see § 5 for discussion). 

Throughout, a ACDM cosmological model is assumed 
with H = 7O/170 km s^Mpc -1 , fi m = 0.23, Q b = 0.039 
and Oa = 1 — ^m- Due to the small redshift range of 
the observed clusters, we neglect redshift dependencies 
where justified. 

2. OBSERVATIONS 

We use the compilation of Rephaeli et al. (2008) for 
the RXTE, BeppoSAX and INTEGRAL observations, 
the results of Ajello et al. (2009) for the Swift/BAT ob- 
servations and the results of Wik et al. (2009) for the 
Suzaku observations. We consider 6 clusters observed 
with RXTE (in the 20 - 80keV band), 4 clusters ob- 
served with BeppoSAX (20 — 80keV), 1 cluster observed 
with INTEGRAL (44-107keV), 8 clusters observed with 
Swift/BAT (50 - 100 keV) and 1 cluster observed with 
Suzaku (12 — 70kcV), where the Coma cluster has been 
observed with all five instruments, A2319 with three and 
A2256 with two. Our sample includes therefore a total of 
13 different clusters, for which HXR detection or upper 
limits are available. Most of the clusters show clear signs 
of a recent merger (this may bias the inferred model pa- 
rameter values, see discussion in § 5). A list of all the 
clusters in our sample, including their relevant proper- 
ties, is given in table 1. 

Different instruments are sensitive over different X-ray 
energy ranges. When comparing the measurements of 
different instruments we therefore compare the flux per 
logarithmic energy interval of X-ray photons, i.e. the flux 
divided by A = log(z/ max /f min ) where hv max , min are the 
upper and lower bounds respectively of the instrument's 
energy band. For the relevant instruments, A is within 
the range 0.7 — 1.8. Our choice is motivated by the fact 
that a photon spectral index of 2, i.e. dn 1 jde 1 oc e~ 2 , 
is expected in many models (see § 3 & § 4) and is also 
consistent with observations (note, however, that there 
are large uncertainties in the observational determination 
of the spectral index, e.g. Rephaeli et al. 2008). For such 
a spectrum, the flux per logarithmic energy interval is 
independent of energy. 

3. A SIMPLE MODEL FOR THE HXR EMISSION 

Let us first briefly describe the main assumptions, and 
the main relevant results, of the model discussed by 
Kushnir & Waxman (2009) for the nonthermal emission 
produced by cluster accretion shocks. In this model, it 
is assumed that matter is accreted onto a cluster of mass 
M at a rate M = /i n st^20o/te , where tn is the (in- 
stantaneous) Hubble time, M200 is the mass contained 
within a radius r2oo, within which the mean density is 
200 times the critical density p C rit, and /; ns t is a dimcn- 
sionless parameter of order unity, reflecting the temporal 
fluctuations of M j (Mzoo/tu)- As discussed in Kushnir & 
Waxman (2009), 3D numerical simulations indicate that 
the average value of /i nst is ~ 0.5. The accreted gas is as- 
sumed to be shocked to the cluster's virial temperature 
T (see discussion in § 5). Since the accretion shock is 
strong and collisionless, it is assumed that it produces a 
nonthermal population of relativistic electrons, with an 
energy spectrum (e.g., Blandford & Eichler 1987) 

dn/de oc e~ 2 . (1) 
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Table 1. The properties of clusters included in our sample 



Cluster name 


T [keV] 





2 


»*200 [h 7 Q M P C 1 


r c [h n l kpc] 


L X [10 45 h^crgs- 1 ] 2 


HXR Flux [10- 12 crgcm- 2 s" 1 ] 


Coma" 


o OO+0.34 
— 0.34 


0.654+ /"? 

— 0.021 


0.0232 


2.30 


246 


1.14 


21 ± 6 R 
15 ± 5 B 
18 ± ll 1 
1.7±0.17 Sw 
< 10 Sz 


A2319 a 


8.8 ± 0.5 


n koi +0 013 


0.0564 


2.26 


204 


2.44 


14 ± 3 R 














< 23 B 
< 0.67 Sw 


A2256 a 


6.6 ±0.4 


0.914lg;gJ* 


0.0601 


2.40 


419 


1.16 


4.6±2.4 R 














8.9+ 4 ° B 


A2163° 


13.29 ±0.64 


n vnc+0.03 

0.796I 028 


0.201 


3.21 


371 


6.62 


11+17 R 
ll -9 


A3667 a 


7 ±0.6 


0.541 ± 0.008 


0.056 


1.92 


199 


1.24 


<4 R 


1ES0657-55.8 6 


17.4 ±2.5 


0.62 ± 0.007 


0.296 


3.22 


257 


7.14 


5±3 R 


A2199 a 


4.1 ±0.08 


655+ ' 019 

U.DOO_ () Q21 


0.0302 


1.62 


99.3 


0.404 


9.8 ±4 B 


A3266 a 


8 ±0.5 


U- '»D_o .019 


0.0594 


2.46 


403 


1.22 


< 0.57 Sw 


A3571 a 


6.9 ± 0.2 


0.613 ±0.01 


0.0397 


2.04 


129 


1.04 


1.4±0.5 Sw 


A2029 a 


9.1 ± 1 


0.582 ± 0.004 


0.0767 


2.29 


59.3 


2.62 


< 1.27 Sw 


A2142 a 




0.591 ±0.006 


0.0899 


2.36 


110 


3.36 


< 1.50 Sw 


Triangulum™ 


9.6 ± 0.6 


0.61 ±0.01 


0.051 


2.39 


199 


1.96 


< 0.65 Sw 


Ophiuchus" 


10.26 ±0.32 


74? +0.035 
u - -0.032 


0.028 


2.73 


199 


1.95 


< 2.80 Sw 



1 Error bars given when relevant to our analysis. 
2 Bolometric luminosity 

a Data taken from Rciprich & Bohringcr (2002) 
b Data taken from Tucker ct al. (1998) 
B BcppoSAX measurements 
INTEGRAL measurements 
R RXTE measurements 
w Swift/BAT measurements 
,z Suzaku measurements 



The fraction of the post shock thermal energy density 
carried by relativistic electrons is denoted by rj e . 

The accelerated electrons lose energy by IC scatter- 
ing of CMB photons (which dominates over synchrotron 
emission at the accretion shock). At sufficiently high en- 
ergies, where the electron cooling time is short compared 
to the time scale for cluster evolution, the resulting IC 
luminosity per logarithmic photon energy interval is sim- 
ply given by 



Al 



200 



^£lC,shock _ 1 3 rj e f b T 

2 2 A e [irriptH ' 



(2) 



where fb = 0(,/fi m , \im p is the average mass of shocked 
plasma particles, and A e ~ 20 is the number of loga- 
rithmic energy intervals in the energy spectrum of the 
relativistic electrons. The frequency at which an elec- 
tron emits most of its IC power is given by v = vo"/ 2 , 
where v$ = 'STcuB/h, Tcmb is the CMB temperature 
and 7 is the Lorentz factor of the electron. The Lorentz 
factor of electrons emitting HXR photons of energy £hxr 



is 



Thxr « 5.3 x 103 (^) 
and their cooling time is 

-1/2 



1/2 



0.44 



20 keV 



Gyr. 



(3) 



(4) 



Thus, the cooling time of electrons emitting at the HXR 
band is short compared to the clusters' evolution time, 



and the IC luminosity at the HXR band is well approxi- 
mated by eq. (2). 

In order to determine the cluster's HXR luminosity 
(and surface brightness) as function of its temperature 
T, a relation between T and M 2 oo (and r2oo) should be 
used. The density profile of the X-ray emitting ICM 
is generally well described by a "/3- model" (Cavaliere & 
Fusco-Fcmiano 1976; Gorenstein et al. 1978; Jones & For- 
man 1984), 



PgasM oc 1 + — 



-(3/2)/3 



(5) 



where r c is the X-ray core radius. Assuming the ICM 
plasma to be isothermal and in hydrostatic equilibrium 
gives 



7"200 ' 



/800tt \ 1/2 / 3/3T 

'3.1^ 1/2 r i 1/2 ^ 70 1 Mpc, 
/800tt \" 1/2 f_3/3J_\ 3/2 

^3.5- 10 15 /3 3/2 T 3/2 i 



1/2 



200 ' 



(6) 



Here, Ti = T/10 keV and [i ~ 0.59 is the mean molecular 
weight for fully ionized gas with hydrogen mass fraction 
of x = 0.75. We assume that the accretion shock is lo- 
cated at r ~ T200) since spherical collapse models predict 
a cluster virial density < p v ir >— 178p cr ;t for £l m = 1, 
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Q\ = (with weak dependance on the background cos- 
mology for the relevant range 0.3 < Q m < 1). 

In addition to the assumption that the ICM is isother- 
mal and in hydrostatic equilibrium, we have assumed in 
deriving eq. (6) that the density profile is given by eq. (5) 
out to r2oo- This implies p oc r~ 2 at large radii, in con- 
trast with the p oc r~ 3 dependence expected at large r 
(e.g. Navarro et al. 1997). However, a detailed discus- 
sion of the accuracy of cluster mass determination un- 
der these approximations, given in Reiprich & Bohringer 
(2002), shows that eq. (6) may overestimate M200 by no 
more than 20%. 

Since most of the clusters observed in HXR have re- 
cently undergone a merger, a note is in place concern- 
ing the validity of our approximate description of clus- 
ter properties (eqs. 2, 5, 6). X-ray cluster maps (e.g., 
Finoguenov et al. 2005) show that both relaxed and un- 
relaxed clusters are not simple hydrostatic equilibrium 
systems. Particularly, the structure of recently merged 
clusters often deviates from a hydrostatic equilibrium. 
However, since the deviations of the gas profiles from 
the mean hydrostatic equilibrium profiles are at the few 
tens of percent level, we expect our description to be 
approximately valid. Detailed analysis of numerical sim- 
ulations of cluster mergers (e.g. Ricker & Sarazin 2001; 
Poole et al. 2006, 2007) can be used to determine the ef- 
fects of deviations from our simple approximate descrip- 
tion. 

Using eq. (6), the accretions shock luminosity, eq. (2), 
gives 

^C,shock = 1J . 10 44 (/ . nst?7e) _ i/3 3/2 
fb 



0.17 



X, 



5/2 



ergs 



(7) 



where (/i„st'7e)_ 1 = /inst W 10 Assuming that the 
HXR emission originates in a thin shell lying at the shock 
radius, the HXR flux of a cluster at a distance d ~ czj Ho, 
within the energy band [£1,22] and within a disk of an- 
gular radius 9 centered at the cluster center, is given by 

F [£u£2] (0) = 7.5 • 10- 8 ((/inst W)_x /? 1/2 

T x ASa.cc. (0)/i 7 o — T 



0.17 



(8) 



where 



Saec.(#) — 2#2 



(9) 



Here #200 = and we have used (/inst)e, defined 

as the average value of /i ns t over the disk considered, to 
explicitly reflect the possible spatial dependence of the 
accretion mass flux. 

Eq (9) gives an approximate description of the depen- 
dence of F on 0, for the case where the emission takes 
place within a shell of radius r = r2oo and thickness 
w <C ^200- The thickness of the emitting region is ap- 
proximately given by the product of the cooling time 
of the emitting electrons and the velocity of the down- 
stream fluid relative to the shock velocity, Ud- Using 

u d = ^/T/3pm p ~ 7.4 • 10 2 T 1 1/2 kms" 1 , we have 



W ~ U d tcool 



; , 3 . 10 2 Til/2( _™ ) -- kpc _ 



(10) 
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Fig. 1. — -F/Ag acc . as function of T/3 1 / 3 (see eq. 8). A linear fit 
for ln(F/Ag acc .) as function of In^Tp 1 / 3 ) gives a slope of 1.4 ±0.5, 
consistent with the predicted slope of 3/2. For a slope of 3/2, 
the best linear fit is obtained for /j ns tf? e ~ 0.2. In deriving the 
fits, we used average values (over all instruments) for the Coma 
and 712256 fluxes, and disregarded the Swift/BAT measurement, 
for which the existence of a nonthermal component is uncertain 
(Ajello ct al. 2009). The best fit is shown by the solid black line. 
Constant /i ns t5? e lines are shown as dashed magenta lines. 



This confirms that w ^200- For w = 0.1r2oo, eq. (9) is 
accurate to better than ~ 25% for any 0. 

Eq. (8) predicts that F/Ag acc . should scale as 
(Ti/3 1 / 3 ) 3 / 2 , with a normalization that depends on 
{/insOe^e- The model predictions are compared with ob- 
servations in figure 1, where we show the measurements 
and upper- limits of F/Ag acc , for the entire sample, and 
in figure 2, where we show only clusters with HXR detec- 
tions (The vertical error bars reflect only the uncertainty 
in measured fluxes, F; we ignored the errors in the de- 
termination of gacc.) which are significantly smaller). A 
linear fit for ln(.F/A(? acc .) as function of h^T/? 1 / 3 ) gives 
a slope of 1.4 ± 0.5, consistent with the predicted slope 
of 3/2. For a slope of 3/2, the best linear fit is obtained 
for Jinst^e ~ 0.2. The small deviations (up to a fac- 
tor of 2) from a constant fi ns tVe line could result from 
cluster-to-cluster variations of /i ns t , or from variations of 
the accretion flow across individual clusters (note that 
different instruments have different FOV). 

4. OTHER MODELS 

In the preceding section we presented a simple model 
for the HXR emission from galaxy clusters, and have 
shown that it is consistent with the available HXR clus- 
ter data. In this section we consider alternative models, 
in which the HXR emission is predicted to be correlated 
with the cluster thermal emission. In such models, the 
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Fig. 2. — Same as figure 1, including only clusters with HXR 
detections, and average (over all instruments) values for the Coma 
and A2256 clusters. 



HXR emission should be strongly dominated by emission 
from the cluster's core, which dominates the thermal X- 
ray emission. This is in contrast with the predictions of 
the model described in § 3, in which the HXR emission 
is produced at the cluster accretion shock. In the lat- 
ter model, the HXR surface brightness is expected to be 
nearly uniform across the cluster, and enhanced along 
the (apparent) accretion shock ring (see eq. 9). 

A widely discussed model for the HXR emission, which 
predicts the HXR emission to be proportional to the ther- 
mal emission, is a model in which the relativistic elec- 
trons are secondaries produced by inelastic p-p collisions 
between cluster CRs and thermal intra-cluster gas (e.g., 
Dennison 1980). We therefore consider a simple version 
of this model in some detail below. In this simple ver- 
sion, the ratio between the CR energy density and the 
thermal energy density of the gas is constant. Although 
this ratio is predicted to scale with the gas density as 



(Jubelgas et al. 2008; Pfrommer et al. 2008; Kush- 



p -l/3 

nir & Waxman 2009), neglecting the radial dependence 
of this ratio is justified since the emission is dominated 
by the core (see Kushnir & Waxman 2009, for a detailed 
discussion). As explained at the end of this section, our 
conclusion, that the data indicate that the spatial distri- 
butions of the HXR and of the thermal X-ray emission 
are different, is general and applies not only to the sim- 
plified secondary electron model. 

In the secondary electron model, the ICM is assumed 
to contain a population of proton CRs with a power law 
energy distribution e 2 dncn./de = /3 core 3nT/2, where n 
is the ICM number density and /3 C orc is the the ratio 
between the CR energy (per logarithmic particle energy 



interval) and the thermal energy. The HXR luminosity 
per logarithmic frequency interval is proportional in this 
model to the X-ray luminosity (see Kushnir & Waxman 
2009, for details), 



R2 

fl CMB 



1.1-10 /3 coro _ 4 Ti / ^ J^- 

-°CMB ' D 



(11) 



where /? c 



/3corc/10 4 and B is the intra-cluster 



magnetic field. -Bcmb = (S^^cmb) 1 ^ 2 ~ 3.2 /iG is the 
magnetic field, for which the magnetic energy density 
equals the energy density of the CMB. Note, that in 
contrast with the magnetic field at the accretion shock, 
which is much smaller than -Bcmb, the magnetic field at 
the cluster core is expected to be > -Bcmb. 

The flux within a disk of angular radius 6 for a cluster 
at a distance d ~ czj Hq is given (for (3 > 0.5, see Sarazin 
& Bahcall 1977) by 



:5.4- 10~ 15 h% 

■,1/2 



XTj Bx : 45.5Agtherm.(#) 



where 



B_ 2 5 ,B » Scmb a er g 
< 0.1, B«B C mb Pcorc '- 4 cm2s' 



z 



(12) 



Sthc 



= (3/8- 5 



L 



^Coma 
min(6d/r c ,r 2 oo/rc) 



rdr 



(1 + f2) 3/3-V2' 



(13) 



h% L x /3 ■ 10 45 ergs~-, p C ore,-4 



Here Bx,45.5 

A;orc/10~ 4 , and .B_5 = -B/10/iG. In order to obtain an 
analytic relation between the luminosity and the surface 
brightness, we have assumed in deriving equations (12) 
and (13) that the cluster emission extends to infinite ra- 
dius. For (3 values close to 0.5, the derived relation de- 
viates significantly from the one that would be obtained 
assuming that the emission is strongly suppressed be- 
yond r 2 oo- For clusters with /3 close to 0.5 we do not 
use, therefore, eq. (13), but rather the relation obtained 
assuming emission is truncated beyond r2oo- 
In the secondary model, Fj Agthcrm. scales as 
with normalization that depends on 



7,1/2 r 



-B /3 CO rc (or only on /3 core if B < B C mb)- In figure 3 we 

1 /2 

show -F/Agtherm. as function of T x Bx,45.5 (Vertical er- 
ror bars reflect only the uncertainty in measured fluxes, 
F; we ignored the errors in the determination of c/therm., 
which arc significantly smaller). Comparing the detected 
HXR emission with the constant B _2 /3 core (or f3 cole ) lines 
of the figure, we find that Bzl(3 C orc-4 — 10 4 (or /3 C oro,-4 = 
10 3 for B <C -Bcmb) is required in order to account for 
the detected fluxes. This implies that the total energy 
density of CR protons, log(e max /e min ) x j3 COTe m 20/3 core , 
should exceed the thermal energy density of the gas in or- 
der for the emission from secondary electrons to account 
for the observed HXR emission. This strongly disfavors 
the secondary electron model. Moreover, the correlation 
between the radio and thermal X-ray emission of galaxy 
clusters suggests (Kushnir et al. 2009) that the typical 
value of /3coro is — 2 x 10~ 4 and that B 3> Bcmb for 
clusters which host radio halos. 
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Fig. 3. — F/ Ag t ^ CInl as function of T 1 ' Lx,45.5- Dashed ma- 
genta lines show constant values of B~ 2 j3 cove (or /3 C ore, see eq. 12). 



Examining fig. 3, we find that the Swift/BAT detection 
and upper limits on F / Ac/thcrm. are well below the values 
of -F/Agthcrm. inferred for clusters with similar T from 
the detections of other instruments. The systematically 
lower HXR fluxes derived from the Swift measurements 
are naturally explained by the model described in § 3, in 
which the HXR emission is extended and not dominated 
by the cluster core. In this case, the lower Swift fluxes are 
due to its smaller FOV. Figure 1 demonstrates that this 
is a valid explanation (the upper limits on F/hg acc . ob- 
tained by Swift are consistent with the values of F/Ag acc , 
inferred from the detection of other instruments). The 
lower HXR fluxes implied by Swift's observations may 
also be explained by assuming that the clusters observed 
by Swift are intrinsically different (e.g. having signif- 
icantly lower values of /i ns t?7e) than those observed by 
other instruments. Since we have no reason to assume 
that the Swift clusters are systematically different than 
those observed by other instruments (all are prominently 
merging clusters), we conclude that the data disfavor all 
models, in which the spatial distribution of the HXR and 
of the thermal X-ray emission are strongly correlated. 

5. DISCUSSION 

We have presented in § 3 a simple model, that explains 
the HXR emission from galaxy clusters as IC scattering 
of CMB photons by relativistic electrons accelerated in 
the accretion shock surrounding the cluster: The cor- 
relation predicted in this model between the HXR sur- 
face brightness and the cluster temperature is consistent 
with the observations, and the observed HXR luminos- 
ity is consistent with the fraction rj e of shock thermal 
energy deposited in relativistic electrons being rj e ~ 0.1 
(see fig. 2). The implied acceleration efficiency of elec- 



trons is similar to the acceleration efficiency of protons 
in the accretion shocks, which is inferred from the cor- 
relation between the radio flux and the thermal flux of 
galaxy clusters (Kushnir et al. 2009). The nonthermal 
luminosity and surface brightness produced by the accre- 
tion shock are determined in this model by the cluster 
thermal properties, and are given by eqs. (7) and (8). 

Several comments are in place here regarding the es- 
timated value of r) e . HXR observations do not allow 
one to determine rj e directly. Rather, such observations 
constrain directly only the value of the product r] e fi nstl 
where /; ns t is the mass accretion rate measured in units 
of M 2 oo/tH (see eq. (7)). We have found that the ob- 
served HXR fluxes are consistent with rj e fi nst ~ 0.1. 
However, since the sample of clusters for which HXR ob- 
servations are available is not complete, the inferred value 
of f? e /inst ~ 0.1 may be biased, i.e. may differ from its av- 
erage value (over all clusters). In particular, since most 
of the clusters chosen for HXR observations are merging 
systems, in which enhancement of the accretion rate is 
expected (see e.g., Pfrommer et al. 2008), the inferred 
value of ^e/inst is probably biased in the current sam- 
ple towards values higher than average. Determination 
of the average value of /i ns t from numerical simulations 
and using a complete cluster sample, that may be pro- 
duced by future HXR missions (e.g. NuStar, Simbol-X), 
would allow one to estimate the value rj e more accurately. 

Another limitation of the estimate of rj e should be men- 
tioned. We have assumed in our analysis that the ICM 
plasma is isothermal and in hydrostatic equilibrium. De- 
viations from this simple model near the virial radius 
may change our estimates for r/ e . Although such devia- 
tions are only weakly constrained by observations, both 
observational (e.g. Vikhlinin et al. 2005) and theoretical 
(e.g. Roncarelli et al. 2006) analyses indicated that they 
are not large (for example, the accretion shock temper- 
ature is lower than the virial temperature by no more 
than a factor ~ 2). Modifications of the ICM properties 
near the virial radius may be easily incorporated into 
our model. Improved (observational) determination of 
the ICM profile near the virial radius will therefore allow 
one to improve the accuracy of the determination of rj e . 

Our model predicts the HXR surface brightness to be 
nearly uniform across the cluster, and enhanced along 
the (apparent) accretion shock ring (see eq. 9). This is 
in contrast with models, in which the HXR emission is 
strongly correlated with the thermal X-ray emission, that 
is dominated by the cluster's core. We have shown in § 4 
that the low values of HXR flux inferred from Swift's 
observations disfavor models, in which the HXR emis- 
sion is dominated by the cores of the clusters: For an 
extended HXR emission, the low Swift fluxes are natu- 
rally explained as due to the smaller FOV of Swift, while 
for emission dominated by the cores of clusters the flux 
should not depend strongly on the FOV (see last para- 
graph of § 4 and compare figures 1 and 3). Moreover, it 
was shown in § 4 that a widely discussed model for HXR 
emission, in which the relativistic electrons producing 
the radiation are secondaries produced by inelastic p-p 
collisions between cluster CRs and thermal ICM (e.g., 
Dennison 1980), requires the total energy density of CRs 
to exceed the thermal energy density of the ICM in order 
to account for detected HXR fluxes. 

Our model prediction, that the HXR emission is ex- 



7 



tended, may be tested by future HXR missions capable 
of producing high resolution HXR maps of clusters (NuS- 
tar and Simbol-X, e.g., may reach a resolution of tens of 
arcsec). This prediction could also be tested by future 
7-ray observations. Our model predicts that cluster ac- 
cretion shocks produce a 7-ray flux of (see eq. (8)) 



^tfW = 4-7 • 10- 6 ({f^eVeU P 



1/2 



fb 



p 3/2 / £i/,p 

1 V10G 



-r 1 



^0.177 1 VlOGcV/ 
xgacc.W^oPhcm-^- 1 . (14) 

This flux, which would have been marginally detectable 
by EGRET (see detailed discussion for the Coma cluster 
in Kushnir & Waxman 2009), should be easily detectable 
by Fermi, which has a ~ 50 times higher sensitivity and 
which may resolve the cluster (the angular resolution 
of Fermi reaches 0.1° above 10 GcV, see http://www- 
glast.stanford.edu/). 

Imaging Cerenkov telescopes may also detect the pre- 
dicted nonthermal flux. However, it should be noted 
that flux predictions for energies > 1 TeV are uncer- 
tain. ~ 1 TeV photons are expected to be produced by 
the highest energy electrons accelerated in the accretion 
shocks (Loeb & Waxman 2000). Since, however, our es- 
timate of the cutoff energy of the electrons is not robust 
(see Keshet et al. 2003, for details), the > 1 TeV flux 
may fall well below the prediction of eq. (14). Lowering 
the energy threshold of the imaging Cerenkov telescopes 



to ~ 0.1 TeV would be very helpful in this context, since 
the prediction of eq. (14) is more reliable at photon en- 
ergies <C 1 TeV. It is important to mention here that 
measurements of the nonthermal emission in different 
energy bands (e.g. HXR and 7-rays) would allow one 
to constrain the energy distribution of the accelerated 
electrons. 

A comment is appropriate regarding the synchrotron 
emission from the accretion shocks. Since the mag- 
netic field at the accretion shock is expected to be weak, 
~ 0.1 /iG < £>cmb (Waxman & Loeb 2000), the syn- 
chrotron surface brightness produced by the accretion 
shock is negligible compared to that produced by sec- 
ondaries (e.g., Kushnir & Waxman 2009). Detection of 
the accretion shock synchrotron emission is unlikely with 
present-day radio telescopes, but should be possible with 
next-generation telescopes such as the LOFAR and the 
SKA (for detailed discussion, see Keshet et al. 2004). 

We finally note that our estimates for the IC flux can- 
not be used directly for the Soft X-ray band (< IkeV), 
since the cooling time of the emitting electrons may ex- 
ceed the dynamical time of the cluster (see eq. (4)). 
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